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The class of 'Asymmetric Dark Matter' scenarios relies on the existence of a primordial 
particle/anti-particle asymmetry in the dark sector related to the baryon asymmetry as a 
way to address the observed similarity between the baryonic and dark matter energy densi- 
ties today. Focusing on this framework we calculate the evolution of the dark matter relic 
abundance in the presence of particle/anti-particle oscillations. We show how oscillations re- 
open the parameter space of asymmetric dark matter models, in particular in the direction of 
allowing large (WIMP-scale) DM masses. Finally, we constrain the parameter space in this 
framework by applying up-to-date bounds from indirect detection signals on annihilating DM. 

1 Introduction 

In the last years we witnessed a revival of interest in Asymmetric Dark Matter (aDM) models 
(for a recent list of references seeEE) in which dark matter (DM) abundance is determined by an 
initial asymmetry between dark matter particles and anti-particles. The main motivation for this 
approach is to address the observed similarity in the DM abundance (r^DM) and the abundance 
of baryons (J^b); ^b/^dm ~ 1/5 by relating primordial asymmetries in the two sectors'^. 

Like in the case of baryons, if there exists a primordial particle/anti-particle asymmetry in 
the dark sector, annihilations proceed until depletion of the symmetric part of the population. 
The final abundance of DM is therefore controlled by the primordial asymmetry, in a contrast to 
usual freeze-out WIMP scenarios where the relevant quantity is the value of the DM annihilation 
cross section. A possible issue in this approach is that if the asymmetry in the dark sector is 
linked to the baryon number, then the observed DM abundance is explained naturally for rather 
hght, O (5 GeV) DM mass. 

This conclusion changes in the presence of oscillations between DM and anti-DM particles, 
and in we study this process in detail (see alsoS?!). Oscillations replenish the depleted pop- 
ulation of 'targets' so that annihilations can re-couple and deplete further the DM/anti-DM 
abundance. The final DM relic abundance is therefore attained through a more complex history 
than in the standard aDM case and in a closer similarity to the freeze-out one opening up the 
parameter space towards heavier DM masses in a natural way. 

2 Theory motivations 

we assume that the dark matter particle DM is not its antiparticle DM and that there exists 
a primordial asymmetry between the two populations. We then study the evolution of the two 



populations in the presence of oscillations generated by a A(DM) = 2 mass term, 5m. The 
effect of A (DM) = 2 operators is to introduce a mass splitting and mixing between DM and 
DM, which are no longer mass eigenstates. We will use the generic form of a Hamiltonian 

^, I m Sm \ , ^ I A if fermionic DM 

^=[6m m) '^"^=\AV(4M) if bosonic DM 

The value of the 'Majorana' mass A is naturally much smaller than the 'Dirac' mass m 
since A violates a global C/(1)dm symmetry. In our model-independent study, we leave 6m to 
be a free parameter which will be scanned over orders of magnitude in the sub-eV range. A 
natural value for this parameter in the fermionic case is obtained from the dimension-5 operator 
^^j^''^ ■ After electroweak symmetry breaking and taking A at the Planck scale we obtain the 
see-saw value Sm ~ 10~^ eV. This value of 6m turns our to be comparable to the value of the 
expansion rate at the time of WIMP freeze out and it is therefore expected that it will lead to 
interesting cosmological effects, as me demonstrate below. 

3 Oscillation -|- annihilation -|- scattering formalism 

Our aim is to study the evolution in time t of the populations of DM particles ( n"*") and 
their antiparticles DM, (n~), which are subject to the simultaneous processes of annihilations 
DM DM SMSM (with SM being any Standard Model particle), oscillations DM o DM and 
elastic scatterings DM SM — )• DM SM. For definiteness, we assume that particles are initially 
more abundant than antiparticles, i.e. n+ > n~. 

This problem, in which a coherent process such as oscillations is overlapping with incoherent 
processes such as annihilations and scatterings is best described by the density matrix formalism, 
originally developed for the case of neutrino oscillations in the Early Universe^'. In this approach 
a 2 X 2 matrix is defined with diagonal entries correspond to the individual number densities 

and n~ and whose off-diagonal entries express the superposition of quantum states ^ and ~ 
originated by the oscillations. As is customary, we recast the problem in terms of the co-moving 
densities = n^/s, where s is the total entropy density of the Universe, and we follow the 
evolution in terms of the dimensionless variable x = m^yM/T, where m^^M is the DM mass and 
T the temperature. We will therefore work in terms of a comoving number density matrix 

(the curly font for y will indicate in the following the matrix quantity). We will always be 
interested in the epoch of radiation domination, during which the Hubble parameter H{x) = 
^ Svr^f?* (x) / 90 m^^x^^/Mpi = Hm/x"^ and t^^ = 2H{x). Here 5*(a;) and g*s{x) are the effective 
relativistic degrees of freedom. 

The evolution equation for the density matrix y then reads (with ' = ^ j^^^^ ^ 

s{x) 



On the right hand side, the first term accounts for oscillations, the second for annihilations 
and the third for elastic scatterings, with the corresponding operators defined as (for more detail 



see' 



« = ( + "fi + „, V = I '= : I and r.r' = L (4) 




(5m TTT-DM + 

where (civ) is the thermally averaged annihilation cross section and AV represents the effective 
energy shift of DM versus DM induced by the baryon asymmetry of the medium and 7^ is 
the elastic scattering cross section. Both AV and 7^ depend on the scattering cross section of 
DM with ordinary matter and we parametrize them through the ratio with the standard 'weak 
interaction' values. The direct detection experiments constrain this ratio to be < 10~^, and we 
will use two values ^ = and ^ = 10"'^ in the following. The common terms on the diagonal of 
Ti, Dirac mass term rriDM and an effective matter potential V do not affect oscillations. 

In eq. (3), l^q denotes an equilibrium comoving density l^q = (f )^^^ ^x^/^e"^, where 
g is the number of internal degrees of freedom and y is the charge-conjugated matrix of 
i.e. the same quantity as the latter but with the role of particles and antiparticles flipped: 
y = CP-^ • 3^ • CP, where CP = ias = ( J ) . The initial conditions read = Y^{xo) = 

Yeq{,XQ) e^^° and (xq) = Y ''(xo) = 0, at an initial time xq. It is also useful to introduce 
the parameter r/o = ~ > which represents the initial DM - DM asymmetry and is related 
to Co as Co = arcsinh(r/o/(2yeq(a:o)))- 

The interplay of the coherent and incoherent processes (annihilations and scatterings) can 
thus be thoroughly followed by using the full density matrix formalism Q We summarize our 
results by showing the contour lines corresponding to the correct DM abundance in the (m,DM; fo) 
plane in Fig. [Tj By changing the values of 5m (and in general also r/Oi which is not shown here) 
we can open much more of the parameter space, towards larger m-oM and larger ctq. 

As oscillations symmetrize the dark sector allowing for annihilations to resume, the param- 
eter space presented above is subject to the usual constraints on annihilating dark matter. The 
most stringent constraints are set by the indirect DM searches in the present epoch provided 
by the FERMI-LAtHH and by the observation of the Galactic Center halo with the H.E.S.S. 
telescope as well as from considering the effect on the generation of the CMB anisotropics 
at the epoch of recombination (at redshift ~ 1100) and their subsequent evolution down to the 
epoch of reionization'^. We show these constraints superimposed on the available parameter 
space, in Fig. [1} 



4 Summary 

Oscillations arise naturally in aDM models and we have studied the impact of adding oscillations 
between DM and DM particles on the scenario of Asymmetric Dark Matter (see alsoll0). We 
found in particular that a typical WIMP with a mass at the EW scale (~ 1 TeV) having a 
primordial asymmetry of the same order as the baryon asymmetry, naturally gets the correct 
relic abundance if the 5m mass term is in the ~ meV range. This turns out to be a natural 
value for fermonic DM arising from the higher dimensional operator DM^/A where H is the 
Higgs field and A ~ Mpi. 

We have applied the matrix formalism to explore the phenomenologically available parameter 
space, by varying the dark matter mass ?71dm, the annihilation cross section uq, the primordial 
asymmetry in the DM sector r/o and the mass difference 5m, for two discrete choices of the 
parameter C that sets the strength of the elastic scatterings between DM and the plasma. The 

°/n a recent worlJ^ authors derive the density matrix equations for DM oscillations and freeze-out from first 
principles using non equilibrium field theory, and arrive to equations which m practical applications often lead 
to similar results as ours, although in some cases can differ. A paper containing a thorough discussion of these 
issues is in preparation. 
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Figure 1: Summary plots of the parameter space showing the constraints. The dotted lines mark the 
contours along which a correct DM abundance can be obtained. Left: Oscillations and annihilations only 
(S, — 0). Right: Including elastic scatterings = 10~^j. In both panels we assume an initial asymmetry 
Vo — Vb o,nd we show two indicative values of the oscillation parameter 5m. The solid black line at the 
bottom represents the standard case (rj = 0,Sm = 0). The shaded blue regions are excluded by CMB 
constraints, the shaded pink ones are disfavored by gamma ray observations with FERMI and the orange 
ones by observations with H.E.S.S. (see text). The white areas above the solid black line are allowed. 



quantitative results are displayed in Fig [T| together with the most relevant constraints defining 
the still allowed regions. 
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